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means a correspondingly low tensile strength, but if a material could be* 
found with a lower modulus, and a strength not reduced, it would, other 
things being equal, be superior. The search for a material with a low 
dilatation with rise of temperature seems more promising, and it is probably 
on account of its low dilatation that silica brick has proved so useful in 
furnace construction. 
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The present paper is a contribution to the treatment of problems which 
require a solution of the differential equation v 4 ^ = 0. Amongst such 
problems are to be found not only the very slow motions of a viscous fluid in 
two dimensions, but also the flexure of thin flat plates.* 

The prosecution of the investigation has been made possible by the support 
of the Department of Scientific and Industrial Research, which has provided 
financial assistance to enable two of us. to devote the whole of our time to 
the research, and our thanks are offered to the Department for its assistance. 
We also desire to acknowledge the facilities afforded by the Governing Body 
of the Imperial College of Science and Technology in placing a room at our 
disposal in the Department of Aeronautics. 

The method of attack was suggested by the results of an earlier paperf on 
the solution of Laplace's equation, y 2 i/r = 0. It was there found that, for 
any forms of single or double boundary, the solution of problems on the 
torsion of cylinders or the irrotational motion of an inviscid fluid could be 
made to depend on a number of definite integrals. In general, the evalua- 
tion of these integrals involves the use of graphical and mechanical methods, 
and much progress has been made in the direction of simplification of the 
processes. It is hoped that a paper may be presented elsewhere in the near 
future showing the application to a number of engineering problems. 

The solution of the equation V 4 ^ ^ given in the present paper follows 
the earlier work in its generality as to boundary forms and also in the 
general dependence on graphical and mechanical integration for arithmetical 

* Lamb, ' Hydrodynamics,' p. 604. 

+ ' Roy. Soc. Proc.,' A, vol. 95, pp. 457-475 (1919). 
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results. It happens, however, that the sample chosen for illustration, the 
flow of a viscous fluid past a circular cylinder in a parallel- walled channel,, 
admits of the use of purely analytical methods for the greater part of the 
work. It has been found possible to give an approximate expression for the 
stream function -v|r in analytical form, nd „ie stream-lines illustrated in 
fig. 6 were obtained from the formula given on page 409 (equation 57). 

Method of Solution, 

d 2 d 2 
Denoting by y 2 the differential operator k^ + k^ ^ ne solution of the 

equation 

VV = (1) 

is required for certain boundary conditions. In the case of slow motion in a 
viscous fluid such boundary conditions are given in the form 

yfr constant on any solid boundary-^ 

and dyjr/dn zero on any solid boundary J 

These conditions are the mathematical expression of the physical idea that 
there is no relative motion between a solid and the fluid in contact with it.. 
The validity of these boundary conditions has received very strong support 
from recent experiments by Dr. T. E. Stanton, F.R.S.,* in which velocities have 
been measured in pipes at very small distances from the walls, and even in 
the case of violently unsteady motion there is no indication of any failure of 
the equations as now understood to account for observed motion with the 
boundary conditions of no slipping. A critical discussion of earlier experi- 
ments is to be found in the standard treatises dealing with fluid motion.f 
If, following normal usage, we write 

? = V 2 f (3) 

the differential equation (1) becomes 

V 2 ? = 0, (4) 

but the boundary conditions are not in a form suitable for a solution of (4) 
by the methods of the paper dealing with Laplace's equation. From the 
ideas there developed, however, we adopt that which says that a solution of 
(4) may be obtained in terms of the boundary values of J, and that a certain 
distribution of simple sources on the boundary can be found which satisfies 
the conditions of the problem. For reasons connected wholly with the 

* ' Phil. Trans.,' A, vol. 214, p. 199 (1913). 
t Lamb, ' Hydrodynamics, 5 p. 571. 
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mathematical manipulation of the expressions obtained, the simple sources 
have been abandoned in favour of double sources,* and a suitable form for f is 



y 



sin 7 7 



(5) 



where the value of fat any point in the field is regarded as the sum of the 
effects of doublets of strength d% situated on the boundary. The axes of 
the double sources are normal to the boundary ; 7 is the angle between the 
normal to the boundary at the elementary doublet and the line joining this 
doublet to the point at which f is being estimated ; similarly, r is the distance 
between such point and the elementary doublet. 



Application to the Begion inside a Single Boundary. 
In terms of the quantities defined by fig. 1, it is assumed that 



& 



QP Xa 



(6) 



the integral extending round the whole of the boundary. 



Normal at Q 




Fig. 1, 

The scale of the diagram does not enter into the solution of equation (1) in 
any other way than that of similarity, i.e., the values of yjr at homologous 
points are independent of the scale, and hence all quantities representing 
lengths may be regarded as pure numbers expressing those lengths in terms 
of some arbitrary unit. The idea somewhat simplifies the writing of 
formula). 

* Lamb, ' Hydrodynamics,' pp. 55-58. 
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In one of the forms used in the theory of attractions, G-reen's theorem 
states that 

27n|r M = (T & log PM dajp dy P + f ^ d0 M - j ^ log FM <&* (7) 

where the line integrals in question are to be taken round the boundary. In 
this equation it may be noted that the two line integrals are immediately 
calculable from the boundary conditions appropriate to viscous fluid motion,, 
and that th* only undetermined quantity when M is on the boundary is £p. 
Equation (7) is then an integral equation for the finding of £p from the known 
boundary conditions. 

Using the value of £p given by (6) in (7) leads to a modified result in 
terms of boundary doublets d^q \ in the process of substitution, however, a 
change in the order of integration has been made to obtain 

27n|r M = [-[[[ ^jp log PM dstk dy v \dx q + f ^f d0 ¥M - j gj log FM ds ¥ . 

(8) 

The value of I Kp" log PM dx P dy v depends only on the geometry of the 

boundary and the positions of the points Q and M, and may therefore be deter- 
mined in any case without prior knowledge of the distribution of doublets 
appropriate to the physical problem. By a second use of Green's theorem, 
this evaluation is made to depend on single integrals, for 

2tt . QM . log QM . sin y QM = 2 f [ S1I ^ P . log PM . dcc P dy P 

-f J QF . log QF . sin 7qf d0 ¥M ~- j ^- (QF . log QF , sin y qT ) . log FM . ds F 

(9) 
is generally true for the conditions for which Green's theorem is valid. 

Combining equations (8) and (9) so as to eliminate the double integral of 
the latter leads to 

2 7n|r M = U 7T . QM . log QM . sin 7qm — \ QF . log QF . sin y qF d0 mi 

3 



+ i f 9^ (Q F • log Q F . si » 7qf) log FM . ds F \ d Xq 



4- [ f F d0 ¥M -~ f -^f • log FM ds v . (10) 
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It is convenient for further analysis that the differentiation normal to the 
boundary be performed generally (see fig. 2). 

^Yqf _ sin 7 FQ 9 . QF __ 



on F 



dyQv 
ds F 



171 > 



QF 

COSJTfq 

^ QF ' 



3% 

a. qf 

ds F 



■cos7 FQ , 



= — sm7 FQ . 



(11) 



Note that 7 FQ -=f=. y qF . 
3 



dn F 



(QF . log QF . sin 7 QF ) 



djov 



3.QF 



= (l + logQF)QF.sin 7QF ^-logQF . oos 7qp ^£.. (12) 

OS F ' ' CS F 




Normal al F 



Fig. 2. 



Using the values found in (12) to transform (10) gives to the latter the form 
2tt^ m = -Ttt . QM . log QM . sin 7 QM — | QF . log QF . sin 7 QF d0 FM 



+ 1 



i a 



+ log QF) . QF . sin 7 QF . log FM . dy QF 



|- log QF . cos 7 QF . log FM . d F QF I d Xq 






f F ^FM- 



H 



m 



- . log FM . ds F . 



(13) 



F 



If M be brought on to the boundary (13) becomes a modified integral 
equation in which %q is the only unknown. Once % Q has been found to 
satisfy the boundary conditions, equation (13) can be used to find the value 
of ^|r M at any point in the region considered. 
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The Determination of xo,- — It is sufficient to direct attention in the first place to one 
term only of (13) which may be denoted by 



^ M = | QM . log QM . sin y QM dx 



(14) 



in order to see the method of solution for xo,- The point M being restricted to lie on the 
boundary, ^'m may be regarded as a known function of position on that boundary : such 
position will be defined by a length of perimeter " s " measured from some fixed point. 
The first operation is to differentiate in (14) with respect to movement of M along the 
boundary 



31 



d 



dsM 



(QM . log QM . sin y QF ) c^ Q . 



(15) 



Using expressions similar to those in (11) equation. (15) is readily transformed into 



a* 

3«M 



(log QM . cos /3 - sin y M q sin y M q) d^l. 



(16) 



As Q passes through M, fig. 3, the angle /3 passes through zero and cos/3 is sensibly equal 




Fig. 3. 



to unity for small values of QM. On the other hand, log QM is then infinite. The 
expression siny M Q . siny QM is equal to - 1 for small values of QM and the existence of 
this term does not affect the method of solution. 
Now assume that 

Xq = J log QE d&E (17) 

where E E is a function depending on the position of a point E on the boundary. Sub- 
stituting for xo, from (17) and changing the order of integration, equation (16) becomes 






{ 



(log QM . cos /3 - sin y MQ . sin y QM ) d q log QE l cte E , 



(18) 



The internal integral of (18) is discontinuous, for as E passes through M there is an 
abrupt change in the integral from -tt 2 /2 to +tt 2 /2, no matter what the shape of the 
boundary. To see this, consider the integral 



I=JlogQM.c£ Q logQE 
in the case where M and E are fixed, fig. 4. 



(19) 
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In general, for any choice of the point E, two positions Q and Q t can be found such 
that Q X E = QE. Taking two elements of (19) together leads to 

/» 

I - J log |~ <k log QE, (20) 

where the limits are from QE = to QE =• a maximum value. 

As the distance ME is decreased it will be seen that the value of log QM/Q X M tends to 
zero for all pairs of elements at a distance from M large compared with ME. In the 
limit, when E is indefinitely near M, the sensible elements will all be included in a range 
of perimeter for which the angle EMQ is small and the points M, E and Q may be taken 
as lying on the circle of curvature at M. 

If Q, lie between M and E ^_ L (21) 

%1 Q 1 M = ff(:-QEJ v J 

If Q, lie to the left of M * >, k (22) 

u QjM = QE - ae 1 x J 

Owing to change of form of expression, the integral I must be considered in two parts 

QE from to ae 
and QE from ae to aexoc, 

and the expression becomes 



ae cu X oo 



I = f qi^SI^H m loe ?m^ d « qK (23) 

"ae 

Write x for QE/ae in the first integral and \\x for QE/ac in the second. The two 
integrals are found to have the same value and are independent of the radius of 
curvature " a" Equation (23) becomes 



1 = 2 



I log \±^ dr, (24) 

X \—3& 







I is a definite integral of the Gamma function series and its' value is 7r 2 /2. Had E been 
on the other side of M the sign of I would have been changed. 

Owing to the discontinuity in I as shown above, the form of (18) can be changed by a 
partial integration with respect to variation of E ; the result lends itself to an expansion 
for E. The form is 

-7T 2 /2 

l2 = 7T 2 Em + I E dm < \ (log QM cos (3 - sin yqu . sin yqm) dq log QE > , (25) 

ir*/2 

and Em is found by successive substitution. 

In order to shorten the notation when returning to equation (13) we write 

ot m for 2tt^ m — J \jr ¥ dO FM + T^log FM . ds T (26) 

and 

Y MQ for it . QM . log QM . sin y qm — § J QF . log QF . sin y qF . d F # PM 

+ -J | (1 + log QF) QF . sin 7qf . log FM . d F y QF 

-J J log QF . cos 7qf . log FM . d F . QF, (27) 
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so that equation (13) becomes 

f*M = | Y M Q^%Q, (28) 

and combining this with (17) 

w M = J { J Ymq . d q log QE} ^3 E , (29) 



MQ 



dm _ f r ray 

ds M - J i J 3s M • "* " 6 ^ / 
and by the use of a theorem similar to (25) 

3gT ___ 7T 3 

35 M ~~ 2 



. cZ Q log QE ) dH E> (30) 



W 4- 



>— * •» t 1 ^ x MQ 



E< M %^^- lo gQ E f- ( 31 ) 



7T S /4 t_ J ^ M 



Eearranging the terms 



tt 3 ^ dm , j— 4 /* , 

77 #m = — 51 r I fl E <^e 



^-X log QE i , (32) 



and by successive substitution 3m may be found. 

The series for 3 M must of course be convergent if it is to be useful, but 
no proof of convergency has been attempted. 

This completes the general account of the analysis. 

Scope of the Application of the preceding Solution of the Equation y 4 ^ = 0. 

The physical problems which can be solved by the method outlined up to 
the present point appear to be those of the flexure of a thin flat plate 
internal to a given boundary. This boundary may be of any shape in three 
dimensions except for the limitation that the co-ordinate normal to the mean 
plane of the boundary must be small. The slope of the plate at the boundary 
is also arbitrary within the same limits. By means of the theorems relating 
to particular integrals, the solution of the flexural problems may be extended 
to the further stage of an arbitrary loading over the enclosed area. 

If an attempt be made to apply the analysis to the part of infinite space 
external to the boundary, difficulties are encountered and the analogous 
problems in the slow motion of a body through a viscous fluid are included 
in this class. It has already been pointed out by Sir G. G-. Stokes and 
Prof. Lamb* that no solution appears to be possible for the steady motion of 
a cylinder through an infinite expanse of fluid stationary at infinity. It is 
not difficult so to modify (8) that the condition i|r = is satisfied at infinity 
and to examine the equation for the value of %. Failure — in the case of a 
circular cylinder — occurs by numerical equality, but with incompatible signs 
in one or other of the quantities ijr or d^rjdn on the boundary. 

It appeared from considerations of the above difficulties that progress with 

* Lamb, * Hydrodynamics,' pp. 603 and 604. 

VOL. C. — A. 2 D 
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the fluid motion problems would require an extension of the analysis to 
double boundaries, and the method previously found effective for the solution 
of v 2 ^ = in an earlier paper* was again tried and has led to the solution 
given below. The increased labour involved in dealing with double boun- 
daries instead of single ones is considerable. 

Slow Motion of a Viscous Fluid through a Parallel-walled Channel and past a 

Circular Cylinder. 

It has not been found to be necessary to restrict the outer boundary to finite 
distances from the cylinder at all points, and in the direction of flow infinite 
length is retained. 

The details of the slow motion in a parallel channel are well known and 
may be found in standard treatises,*} - and the molecular rotation found tp be 
proportional to the distance from the centre line of the channel. Since (6) 
is only valid for a molecular rotation which vanishes at infinity, this uniform 
streaming has been eliminated from the problem as an initial step and the 
difference due to the cylinder separately estimated. As solutions of the 
equation v 4 ^ = are additive, the two problems are very simply related. 

The procedure, then, is : — 

(1) Ignoring the parallel walls of the channel, find a solution of y 4 ^ = 
which satisfies the required boundary conditions at the surface of the 
cylinder. In the case of the circular cylinder this is readily obtained by 
expressing fin spherical harmonics, but in any other case the method indicated 
in equations (26)-(32) is available. In general it appears that i/r will tend 
to become infinite at an infinite distance from the cylinder. 

(2) Find the values of yfr and <fy/dn from this solution at the position of 
the outer boundary and ignoring the inner boundary, apply equations 
(26)-(32) to find a solution of vV = which has equal and opposite values 
f or i/r and dty/dn at the outer boundary. 

(3) Combine the solutions of (1) and (2) by addition so as to give ^ and 
d^rjdn on the outer boundary. In general the conditions on the cylinder will 
not then be satisfied (in our example a very close approximation is obtained 
at this stage), and the difference between the required and calculated values 
forms a new starting point for a repetition of the process. 

The process lacks the precision of equations (26)-(32), but appears to be 
quite satisfactory in application ; our present example supports the earlier 
one based on vV = in indicating convergency. Connection between the 
two is not wholly lacking since the molecular rotation satisfies the equation 

* ' Roy. Soc. Proc./ A, vol. 95, p. 464 (1919). 
t Lamb, ' Hydrodynamics,' p. 576. 
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V 2 £~ 0, and it is the strength of the f doublets which is found by the 
present analysis. 

The remainder of this paper is devoted to a detailed consideration of the 
solution of a particular problem. 

Initial Conditions. 

A cylinder of unit radius is fixed relative to the parallel walls of a channel 
so that its centre is five units away from each wall. Fluid is forced through 
the channel in quantity sufficient to give unit velocity to the fluid in the 
centre. of the channel at infinity. Absence of relative motion between the 
fluid and solid boundaries is taken as a condition of viscous fluid motion. 
Incompressibility of the fluid is assumed. 

The differential equation of motion to be satisfied is 

vVi = °* (33) 

Along the cylinder 

f 1 = and ibh = 0. (34) 

on 

Whilst along the outer boundaries 

yfn = ± a and ^ = 0, (35) 

on x ' 

where a is determined from the width of the channel and the hypothesis that 
at infinity, where the flow is laminar, the velocity at the centre is unity. 

Corresponding with the special conditions of the problem, we have at 
infinity 

f 2 = _y|l _£/!)*! ( 3 6) 



where 25 is the distance between the two walls, i.e., is equal to 10 (see fig. 5). 
Writing a|t + ^ 2 for -fa, equation (33) is replaced by 

VV = 0, (37) 

yfr = and dyfr/dn = on each of the walls of the channel. 
On the cylinder it is easily found that 



, • A 1 1 sin 2 6 X % 
i|r = sin 6 ( 1 — 



3b 2 J 



on \ b* 



y. (38) 



j 



For the cylinder the normal into the fluid has been chosen in fixing the 
sign of n. 8 is measured from the centre of the cylinder and a line midway 
between the walls. 



2 D 2 



404 Prof. L. Bairstow, Misses R M. Cave and E. D. Lang. 



Ignoring the conditions at the outer boundary, a solution of (37) compatible 
with (38) is obtained by expressing f in terms of spherical harmonics ; it being 
assumed that 



, A • n i , B sin 3 6 , E sin , D sin 3 
ty = At sin v . log r + — + — + — 



■ttO 



(39) 



f = (2-pjrsin^.logr + ^l 



From (38) the values of A, B, C and D are found and give to yfr the value 

1 \ sin 9 sin 3 6 L 2\ ,*^ 

over the whole field. The values of ty and d\jr/dn are required for the 
position of the outer boundary. 

The convenient co-ordinates being K and &, the value of ty on the outer 
boundary is 



f ^(h-^\log(K 2 + b 2 ) + (b + ^-\S. 



1 






(42) 



2b) lK 2 + &* (K 2 + b 2 fJ ' S(K 2 -\-b 2 f 

(41) 

In developing the formula (13) the outward normal was used and hence 

dn ob 

for the outer boundary. Hence df/dn is found from (41) without difficulty. 

The next step is to obtain a solution of the equation (13) as applied to the 

parallel walls. Using m M and Y MQ as defined in equations (26) and (27) the 

solution of 

t*m = | Y MQ dx q (43) 

is required when M is on the boundary wall. It is convenient for later use 
to determine the value of w^tt — 2^ m in a form suitable for points in the 
fluid, and two new variables, h and b 2) are introduced and defined in fig. 5. 






b, 




Fio. 5. 
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The value of 

1 f G\jr 



tsr 



7T 



1 (* 1 C H f 

2-iK = iK<Z0 PM H — U-^ logPMefep (44) 

T ~ 7TJ T TTjOTp ° 



is found by taking the two line integrals together, and in general terms 
leads to 

^ - 2 ^ M = (l - gp) |(6 2 + 6) log (K 2 + h + i) - (h + 6) log (K 2 + 61+6)' 

K , . K 

ran J- 



2 (6 1 -6 2 )-2K tan" 1 -i^-tan 



\ 6 2 + 6 6j + bj 

62 + 6 61 + 6 

+ 2 — ~ 2 (45) 

2 6 2 (K 2 + 6 2 "+ 6) 2 6 a (K a + 61 + 6) 

as a close approximation. In calculating numerical values from (41) it 
appeared that certain terms were very small and in order to reduce the length 
of the expressions they have been omitted in this example. 

When M is on the top right boundary — symmetry makes it unnecessary to 
consider the whole boundary in detail — the value of a/t m is to be equal to that 
given by (41). Putting 6i = 10, 6 2 = and 6 = 5 in (45) a value of mjir is 
obtained which by differentiation with respect to K gives Ijir dm/ds required 
in the solution for 3 by the methods outlined in the earlier part of the paper. 
The analytical expressions are 



GT 



IT 



: 0-98 |l5 log T f 9 2 + ? 5 c: +20-2Kftan- 1 ^-tan- 1 ~ 
L °K 2 + 225 \ 5 15, 

10-3 _ 10-2 20 0-3 afi . 

K 2 +25 (K 2 + 25) 2 (K 2 + 25) 3 K 2 + 225 { } 
and 

1 3?_-n.QQ/9/f «-l K -.fo«-l K \_ 20K 



0-98^ 2 tan" 1 ~- tan 



7r3s< M I \ 5 15/ K 2 + 25 



20*6K 40-8K 40 K _ 0'6K u ^ 



(K 2 + 25) 2 (K 2 + 25) 3 (K 2 + 25) 4 (K 2 + 225) 2 ' 

Without evaluating B M an approximation to ^ Q can be obtained by making 
use of the fact that when % Q is linear between two fixed abscissae for Q this 
segment of % Q contributes to the integral jY MQ dQ for a fixed point, M, a 
quantity proportional to the slope of the segment. This enables a polygonal 
form for Xa to be determined such that the integrals for any chosen points, 
M, shall have any required values. In the special case considered the 
abscissae of the angular points of the polygon were taken at k = 0, 5, 10, 15 
and the chosen positions of M at K = 0, 5, 10, 20. Owing to the symmetrical 
form of Xq. the slope of a side in one quadrant determined a side in each of 
£he other three quadrants. The contributions to the integral due to each 
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group of four sides were calculated for the chosen positions of M on the 
hypothesis that the sides in one quadrant had unit slope. Four simultaneous 
equations then determined the multiples of these contributions — and therefore 
the slopes of the sides — necessary to give sr M the required values at the four 
positions of M selected. In the special case considered the contributions for 
unit slope were calculated analytically, but in general they would have to be 
ascertained graphically. In each quadrant the polygonal form determined as 
above approximated to a straight line. In order to avoid discontinuity the 
straight line was replaced by an arc of a parabola from k = to k = 2*5, and 
then by the tangent to this parabola. For positions of M other than those 
selected the values of m M corresponding to the ^ Q thus determined were fairly 
good approximations to those required, but in order to get a closer agreement 
a small additional term of the type 2a/a 2 + h 2 was introduced and the 
constants were suitably adjusted ; the form of the modifications from the 
straight line was determined by considerations of symmetry and ease of 
analytical manipulation. 

A value of % Q which satisfies (43) very closely is given by 
%Q = * 010 0-Q5937 p ? ic between and 2*5 

iC + £& 

and v = ^ 6 ' 15 -Q-05937 (5&-6'25), k between 2-5 and 0. (48) 

The degree of closeness of the approximation can be seen from Table I 
which shows the value of — m/ir as calculated from (48) and (45) in com- 
parison with the value required. The form chosen is integrable generally for 
wjir in the fluid. The approximation could be further elaborated if required. 

The finding of ^ Q as given by (48) is the essential part of the solution of 
the second stage of the problem of the flow of fluid past a circular cylinder in 
a parallel-walled channel. Corresponding with (48) is a somewhat lengthy 
expression : — - 



^ = 2-615 

7T 



- K 2 -f&! + 5 K 2 + &! + 15 4Ei (h + 5) 45 2 (5 2 + 5) 
log - .. _|_ - 



, 0-05937 
+ ~~6~~ 



K 8 + fc+5 K 2 + 5 2 + 15 K 3 + &i + 5 K 3 + & 3 +5 

2 2 



^+^.f(2-5 + K +9& 2 2 )log(2-5 + Ef + h 2 ) 



2 



-(2-5 + K + 96 1 3 )log(2'5 + K + &i 2 ) 



2 2 



+ (2-5 + K -3 . 6 2 + 10)log(2-5 + K + Z> 2 + 10) 

2 2 2 



(2-5 + K -3 . 6i + 10 ) log (2-5 + K + &1 + 10) 
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+ 2 ' 5 ,. K i (2-5-K + 95 2 2 ) log (2-5-K + &/) 



-(2-5-K + 9V) log (2-5-K -h 2 ) 

2 2 



+ (2~5-~K -3 . & 2 +10 )log(2-5-K +6 a +10 ) 

2 2 2 2 """n 

-(23T^K -3 . & 1 + 10)log(2 T o'^K + h + 10) I 



2 



+ 6»(5V-3.2-5 + K)tan 



, 2-5 + K 



Ji(5&i 3 -3 . 2-5 + K ) tan" 1 2 ' D + K 



2 



2*5 4- K 
— (& 3 + 10)(J 2 -f 10 —3 . 2-5 + K) tan"" 1 ^ , 1Q 



+ (&!+ 10) (6i + 10 -3 . 2'5 + K)tan- 1 -^ ,-T- 

#1 + lu 



+ & 2 (5J 2 2 -3 . 2-5"- K ) tan" 1 2 ' 5 , K 



2 



5 X (55/- 3. 2-5-K) tan" 1 



2-5-K 



2 



(& 2 +10)(& 2 +10 -3 . 2-5-K) tan" 1 , =- 





2-5-K 



2+10 

2 2 



+ (ii + 10)(6 1 + 10 -3 . 2-5-K)tan" 1 



-200(6i-5 2 ) 



2-5-K 

"61+IO 

(49) 



Table I. — Outer Boundary. 



K. 


*_ w 


From (48) 


From (41) 


From (45) and (48) 


— OTM /«*. 


— 5Tm/t. 


3"^/3»M. 


B^/8^m. 





12-30 


12 -30 


5*06 


5-07 


5 


8*39 


8-39 


— . 


4*80 


10 


4-49 


4'53 


5-12 


5 *12 j 


15 


2*62 


2 65 


— 


5*60 


20 


1-75 


1*68 


6-05 


6 *04 


30 


0-83 


0*82 


6 73 


6 *73 


40 


-— - 


— 


7*27 


— 


50 


0*33 


0-33 


7 70 


_ 



The next step is the calculation of ^ M and dty/dn M from (45) and (48) when 
the point M lies on the inner circular boundary. Points being defined as 
before by the angle 0, measured at the centre of the circle, the result of 
the calculation is given in Table II. Since it is of interest to examine the 
effect of molecular rotation, the elements of *\Jr M and d\frfdn M which depend 
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on ^ and those which are independent of it are separately tabulated as well 
as the sum. 





Table II.- 


—•Inner Boundary (d 


ae to Sources on Outer Boundary). 


0°. 


Froin (48) 


From (45) 


^M. 


From (48) 


From (45) 


d^ld^M. 


— STm/t. 


2\|/ — stm/V. 


— 1/7T 30/3»M. 


8/8w M (2j(/-o/7r). 























10 


0-73 


1-57 


0*42 


0-693 


1-576 


0*441 


20 


1-40 


3-09 


0-85 


1-358 


3-104 


0-873 


30 


2-05 


4*52 


1-25 


1-994 


4-500 


1-253 


40 


2-67 


5-81 


1-57 


2-570 


5-790 


1-610 


50 


3-17 


6-92 


1-79 


3*048 


6-832 


1-892 


60 


3*59 


7*82 


2-11 


3-460 


7-780 


2-160 


70 


3'90 


8-48 


2-29 


3 784 


8-489 


2*352 


80 


4-09 


8-88 


2-40 


3-982 


8-830 


2 *424 


90 


4-16 


9'02 


2 '43 


4-035 


9-007 


2-486 
- 



Within the limits of accuracy of working, the above figures are repre- 
sented by 

^ M = 244 sin + 0'Ql sin 30" 

|^ = 2-48 sin + O02 sin 30 ^ ^ ' 



and 



i.e., the values are very nearly proportional to sin 0, with the constants of 
proportionality almost the same. 

Very little further work completes the calculation to the full degree of 
accuracy hitherto attempted. Two values of yfr have been found which 
satisfy the differential equation \7 4 ^ = 0, but with different, though related, 
boundary conditions. As in (38) and (40), define one of these by yjr, and as in 
(44) and (45), define the second by t/r M . The first represents the solution 
which satisfies the inner boundary conditions and ignores the outer boundary, 
whilst the second has the same values as the first at the outer boundary, but 
has not involved direct consideration of the inner boundary. 

On the channel walls therefore 



dn M 



en 



^M = ^. 



(51) 



On the inner boundary, close approximations are 

f = 0-990 sin + 0*003 sin 36h 

= 2-440 sin + 0-010 sin 30 J 



8^ 



(52) 



dn 



0-970 sin + O<OlO sin 30 



-N 



9±_9- 



>■ 



(53) 



dn 



2-480 sin + 0-020 sin 30 



M 



y 
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The function 7 (^ M — ^) can, by choice of 7, be made to satisfy the 

boundary conditions on the cylinder approximately, whilst those on the 

channel walls are completely satisfied. It is, however, easy to carry the 

accuracy further by the use of spherical harmonic functions which vanish at 

infinity and are of the order of accuracy of working on the outer boundary. 

It is assumed, therefore, that 

r / 1 f\. sm # 1 a sin 30 / K /fX 

Y3 = 7(^m — ^) + ^ + p — 3— (o4) 

is to be made to approximate to the boundary values as closely as possible by 
choice of 7, «, and /3. The algebraic work is simple and 

f 3 = 0-662 (^ M -i|r) + 0-030 EEJ?_ 0-001 ^J^, (55) 

gives boundary values on the cylinder of 

^ 3 = 0*990 sin + 0*004 sin 30 " 

and ^p = 0-970 sin 4- 0-010 sin 30 f ' ( ' 

on, 

J 

and, by comparison with those postulated for yfr (see (52) and (53)), the 
accuracy is seen to be complete to the order attempted in the previous Tables. 
On the outer boundary the effect of the added spherical harmonic terms is 
small and does not exceed 0*006 in the value of ^ at any point. 

To complete the solution, it is necessary to add that corresponding with 
steady streaming of the fluid through the unobstructed channel, as given by 
(36), i.e., the final stream function is 

(Qn sin #_sin 30 ( A 2^ 



0-662 (y!r M -~l'96r sin log r- 0*990 ^-E^ (l~ 
{ r lOOr \ or 



+ '030 ^-O-OOl gj^-r sin (l-^ S *f g \ (57) 
r r 3 \ 75 / v 7 

and yfr M is defined by equations (45) and (48). 

From (57) were calculated values of the stream function, from which fig. 6 
was prepared. 

The value of the stream function given by (57) applies to the case of the 
flow of fluid through a channel past a circular cylinder. There is a slightly 
different flow when the cylinder is moved at uniform velocity through 
stationary fluid, and it is perhaps worth while to deduce an approximate 
expression for -fy in the latter case. Strictly calculated, the new problem 
calls for a solution of the differential equation (37) with boundary values 

\|r = sin 0^ 

— -*- = sin 6 ! 



a 



n 
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instead of (38). It is, however, evident that the differences between the two 
sets of conditions are very small, being less than 2 per cent, for ty and 4 per 

///////////////////////////////////////////////////////////////////////////////A . 

"3r3 




Fig. 6.- — Stream lines (^). 

cent, for dyfr/dn. It might then be expected that many of the succeeding 
expressions would differ by similar small amounts. This is seen to be true 
for the equivalent of equation (40), which is 

sin# 



-*fr = 2r sinolog r 4- 



r 



(40a) 



On the cylinder. walls r is unity, whilst the least value on the channel walls 
is 5. The same reasoning leads to the hypothesis that ^ M is sensibly the 
same as that given by (50), and in producing equation (58) the small terms 
in sin 3 in equation (50) were ignored. Following the method indicated in 
equations (54) to (57) a new value for the stream function is obtained as 



f 4 = 0*676 |^m ~l*96r sin <9 log r~ 0*990 



sin0 



+ 



0-020 sin 6 



T 



(58) 



and the value of all the numerical approximations may be tested a posteriori. 

On the cylinder the values of ^ M and d^rjdn M are given by (50), and putting 

r = 1 in (58) leads to 

fa = sin (94 0-007 sin 30 (59) 
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which is very close to the required value. Similar calculations show that 
the remaining boundary conditions, both on the cylinder and channel walls, 
are equally closely approached. . 

Since the relative motion of the fluid and cylinder when the latter moves 
is very little different from that when the fluid moves, it follows that the 
calculated resistance will he given approximately by the same expression (68). 

Calculation of £ and the Resistance of the Cylinder, 

The value of £ may be obtained from (57) by differentiation, but it is 
simpler in the present instance to determine f M directly from the value of %. 
Before integration the expression for f M is 



& = 0-662 



Bin 7 QM 
QM 



^Xq — 3 V 92 



sin 



r 



0-08 -™T + 0'08r sin 0. (5da) 






The value of o\f dy^ was found analytically and the necessary calcula- 
tions were made for the drawing of fig. 7. 



i.o o»& o.e 




Fig. 7. — Lines of constant molecular rotation (£). 

On the inner boundary it was found that a close approximation to f is 

given by 

Cylinder = --2*32 sin #~ 0*053 sin 30 (60) 

and this value was used in estimating the resistance of the cylinder. 
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The formulae for the pressures were taken from Lamb's ' Hydrodynamics, 
p. 570, and with the fluid incompressible are 

du 



P. 



XX 



•p+2fi 



\ 



Pyy = -p + fy 



dx 

dv 
dy 



y 



(61> 



Pyx — Pxy 



(dv du\ 
\dx dy) 



j 



where p is to be obtained from the relations 



dp __ 

dx 

It is not difficult to show that 






2 - ■= 0-662 



cos 7 QM 



d X 



Q 



3- 



92 



dp 

dy 

COS0 



= +f >* 



62> 



0-08 



cos 301 



3 j>+O'O8reos0, (63> 



QM — r i 

and the evaluation on the inner boundary admits of the expression 

^ = /a (2-32 cos + 0-053 cos 30) (64> 

being used as of sufficient accuracy. 

Since f = 3 2 i/r/3r 2 on the cylinder and the velocities along and normal to 
the surface are zero, it is easy to find the differential coefficients du/dx, etc., 
and hence p xx and p xy . 

The resistance of the cylinder per unit length is then 



K = 



"2tt C2tt 

Pxx cos 6dd + pxy sin cl8. 



(65 











It appears that the integrals contribute equal amounts to the resistance,. 

which is, in numerical form 

E = 2*26 irfi. (66) 

It should here be noted that the cylinder is of unit radius and that the free 
stream velocity in the middle of the channel is unity. The formula for 
resistance is easily generalised, for on the principles of dynamical similarity 

Hd 



P V 2 d 2 



=/(v) 



(67) 



where U is the stream velocity, d the diameter of the cylinder and v the 
kinematic viscosity (pv = fi) and p the density of the fluid. Since E is 
proportional to \x it follows that/(U^/v) must take the form AvjUd where A 
is a constant to be determined from (66). It is then found that 

lid 



pWd? 



7-10 



v 



Vd' 



(68) 
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The limitations imposed on the application of the formula are discussed by 
Lamb and it appears that TJd/v should not greatly exceed 0*2 if (68) is to 
-apply. The resistance coefficient RdfpTJM 2 is 35*5 when.TJd/v = 0*2. 

Observations in a wind channel have not been made at values less than 
Vdfp = 10, and the corresponding resistance coefficient is then 1*5, but is 
varying rapidly in the direction leading to high values at low values of XJdfp. 
It may be noted as of interest that the formula given above (68) would give 
a coefficient of 0*71 at Vd/v = 10, is., about half the observed value. The 
departure from fact is then less than might have been anticipated on mathe- 
matical grounds. 

The resistance of a cylinder in a viscous fluid of infinite extent has been 
worked out by Oseen * using a different equation of motion, but subject to the 
same physical and mathematical limitations as in the present problem. The 
resistance as given by Lamb is 

Hd __ 4tir . v/TJd ,«q. 

P V*<P ~~ 1*309-- log TO/// (W) 

.and when Ud/v = # 2 the value of RdfpU 2 d? is 21*6. It is not possible to use 
the formula for higher values of JJd/p on account of the analytical form, since 
the assumption that Ud/v is small has been used both in forming the equations 
•of motion and also in solving them. It appears to be probable that the latter 
limitation can be removed by methods similar to those of the present paper. 
For the moment it may be noted that the presence of a channel only five 
times the diameter of the cylinder in width has not raised the resistance to 
•double that in a free stream. 

* Lamb, * Hydrodynamics/ pp. 605-607. 



